X 



INPP-UVA-00-02 
April 15, 2000 
hep-ph/0008104 



Modified Wandzura-Wilczek Relation 
with the Nachtmann Variable 

O' 
O 

C\I . Xiaotong Song* 

O 

D . Institute of Nuclear and Particle Physics 

Q 



Department of Physics, University of Virginia 
P.O.Box 400714, Charlottesville, VA 22904-4714 

Abstract 



> 

O 

I 

00 

O 
O 
O 

If one retains /Q'^ terms in the kinematics, the Nachtmann variable ^ 
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Typeset using REVT^ 



The standard Bjorken variable xb = Q^/2P-q is commonly used in the discussion of deep 
inelastic scattering (DIS). However, if one retains M'^/Q'^ terms in the kinematics another 
variable ^ introduced by Nachtmann pi] (cf. Greenberg and Bhaumik p|) 



e = 2x/(l + ^l + 4M2xVQ2 ) (1) 

seems to be more appropriate to describe DIS processes. 

For the deep inelastic polarized lepton-nucleon scattering, the asymmetry depends on 
two spin structure functions gi and g2- The structure function gi{x,Q'^) can be interpreted 
as a charge-square weighted quark helicity distribution in the parton model, and the EMC 
measurement led to a surprising result - the quark spins contribute a very small fraction of 
the spin of the proton - the so-called " spin puzzle " . Since then many theoretical works and 
experimental measurements have been done to solve this " puzzle " Most recently, the 
second spin structure function g2, which includes both twist-2 and twist-3 contributions, has 
been measured with relatively higher precision |p. Using the operator product expansion 
(OPE) approach, one can obtain a relation between gi and g2 

g2ix,Q')=gY'^{x,Q')+g2ix,Q'), (2) 

where 

gY'^{x,Q')^-g,{x,Q')+ f'%^{y,Q'). (3) 

Jx y 

If the twist-3 contribution g2{x,Q^) can be neglected, then Eq.(2) reduces to so called the 
Wandzura-Wilczek relation. ^,0] An interesting question is that how significant is the twist- 
3 contribution in g2. Model predictions (for instance see suggest that the twist-3 

contribution is not small compared to gY^ ■ The earlier data given by E143 Collaboration 
||T0| and most recent data given by E155 Collaboration however, seem to show that 
g2{x, Q^) is close to g^^^x, Q^) and the twist-3 part of g2 is rather small. 



In an earlier unpublished note [TT[], we obtained a modified W-W relation (see Eq.(4) 



below) with respect to the Nachtmann variable ^. The result was used by E143 Collaboration 



1^ . In this brief letter, we present the modified g^^iC) and compare it with most recent 



data. More discussions on kinematical corrections arising from the target mass effect are 
given. 

The modified Wandzura-Wilczek relation is 

g. (^^Q) = -sd^^Q) + ^^(^^Q)lj[K;(^^-^^l-j^^ (4) 

where the kinematic factors Ki^2 and K[ are 

^iKU.^ ) (i + MV/Q')(l + 3MV/Q')' ^ 

The derivation of (4) with (5a-c) is given in the appendix. Several remarks are in order. 

• (i) In the large-Q^ limit, all correction factors Kx{i, Q^), Q^) and K2{^-, Q"^) given 
in (5a), (5b), and (5c) approach unity and (4) becomes 



y \ Jy z j 



(6) 



Considering M'^/Q'^ — > 0, and ^ ^ x, Eq.(6) reduces to the original W-W relation 
Eq.(3). 



(ii) From Eq.(l), one would expect ^min = and ^max = 2/(1 + y'l + AM"^ /Q"^ ) for 
X = — > 1. However, since the true momentum fraction carried by quarks is ^ (if 
quark is massless) rather than x, hence we should take ^max = 1 from the beginning. 
We note that the derivation of Eq.(4) does not depend on the value of imax- 
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(iii) To show the correction effect, we first plot i^i, K2 and K[ as functions of ^ 
for = 3 (GeV)2, in Fig.l. One can see that for ^ = ^ 1, i^i = 1 ^ 0.324, 
K2 = 1 ^ 0.460 and K[ = 1 ^ 0.581. It seems that all correction factors reach 
the maximum at ,^ = 0. However, their combined effect presented in Eq.(4) is not so 
simple. 

(iv) To show some aspects of the correction effect, we assume that gi in the integral 
in Eq. (4) is a constant and define a ratio 

g') = /(i^i, i^2, e, 1, q'), (7a) 

where 



-s^r^^V in) 



g2) Q2 4 ^ K((^,g2)^ 

The ratios R{C,,Q^) for = 3, 5, 10 and 100 (GeV/c)^ as functions of ^ are shown 
in Fig. 2. One can see that the ratio is quite large at low and approaches unity 
when oo. However, the ratio (7a) only provides an incomplete information of 

the kinematical target mass correction to . First, the function gi{yiQ'^) is not a 
constant but function of and secondly, one should take the whole result from Eq.(4), 
not just the second term. 



(v) As pointed out in that the original Wandzura-Wilczek relation was derived 
from the Dirac equation for free massless quarks and no higher twist corrections were 
included. By using the equation of motion with nonzero quark mass and imposing the 



gauge invariance, an improved Wandzura-Wilczek relation is obtained in |T3 

g2{x) = -gi{x) + —9i{y)-^ ^ / — r y , 8 

Jx y M Jx y oy Jx y 

where hxix) is the transverse polarization density and T{y) is related to the multi- 
parton distribution hT{x,x'). The quark mass-dependent term (~ rriq/M) in (8) is 
another twist-2 piece in addition to the usual term gY^{x). The last term in (8) is a 



twist-3 term which is coming from the quark gluon interactions. Assuming the niq/M 
term and twist-3 contribution are small, we expect that a modified version of Eq.(8) 
with the kinematical target mass corrections would be very similar to Eq.(4). 

(vi) Making use of the phenomenologically fitted function to the gi data, the modified 
dV^iO Eq.(4), gY^iO Eq.(6), and the original g^^{x) in Eq.(3) are plotted as 
functions of ^ in Fig.3. The data of g2{x) are taken from E143 |T^ and E155 0. From 
Fig. 3, one can see that the effect of the kinematical target mass corrections is rather 
small relative to the experimental errors. All three g^^ curves seems to be consistent 
with the g2 data. However, more precise data are needed for a significant comparison 
of g^^ and g2- Most recently, two papers |14,15| published on the same topic - target 



mass corrections on the Wandzura and Wilczek relation - which found that target mass 
corrections do not affect the W-W relation (2) if all powers in M^/Q^ are included. 
We do not know, however, if this conclusion holds for relation (8). Anyway, since our 
result (4) holds up to the first power of MVQ^ and target mass corrections have very 
small effect on the W-W relation, hence there is no contradiction between ours and 
theirs. 

(vii) It is easy to verify that by changing variable ^ to a; and defining a{x, Q^) = 



1 -|- AM'^x'^ /Q'^ — 1, Eq.(4) can be rewritten as 

g, {x,Q) = -g,{x, Q ) + ^ ^^^^ q^^^ ^ ^ ^^^^ Q^y^ ^.iv, Q 



(l + a(y))(l + a(y,Q2)/2)2 
This is the result obtained in i|T6|l. 



/ -3a(z)(l + a(z)/2)^7i(^,g^)]. (9) 
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Appendix 

From Eqs. (47) and (48) in Wandzura's paper [0|, for n=2,4,..., we have 



where 



2MV2 



Q2 ^• 

In obtaining Eq.(I.l), g2 has been replaced by , or equivalently, g2 has been neglected. 
In the large-Q^ limit, e ^ 0, (I.l) becomes 
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Prom (1. 2a), one easily obtains 

^(e, Q') = Q') + C %i{y, Q') (large - limit), (I.2b) 

y 

which is the same form as Eq.(3), but with respect to the variable ^. Considering C, x 
in the large-Q^ limit, the relation (1.2b) approaches the original Wandzura-Wilczek relation 
Eq.(3). 

On the other hand, in the large-n hmit, (I.l) becomes 
I di.e{l + -^){l + e{i,Q'))[g,{i,Q') + g^'^{^,Q')]^Q (large - n limit) . (I.3a) 

The main contribution to the integral comes from the large-^ region due to the suppression 
factor It implies that 

^(e, Q') ^ -^l(e, Q') (e - Ua.). (I.3b) 

Considering the large-Q^ limit (1. 2b) and large-n limit (1. 3b), it is naturally to assume 

97"^ Q') = -9i{^, Q') + Q') Z'""' dyfiy, Q') (finite g^), (1.4) 



where K2{^,Q^) and f{y,Q^) are two unknown functions to be determined and they have 
the following behavior in the large-Q^ limit 

K2{C,Q')^i, f{y,Q')^^-^^^. (1.5) 

Contrast K2{C, Q"^) with /(y, Q^), the former is a pure kinematical correction factor and 
does not depend on gi. Substituting (1.4) into (I.l) and neglecting the term, we obtain 

/ d^-Ci^ + -^m + ^e)^ + {l + e)K2{C,Q') [ dyf{y,Q')]^0. (1.6) 
Jo Of n + 6 n + 1 Je, 

Unlike the derivation of (1.2b) from (1.2a), we have to use one equation, (1.6), to determine 
two unknown functions. Since K2 is a pure kinematical correction factor satisfies the large-Q^ 
behavior (1.5), we may choose 

K2{C, Q') ^ (1 + 6(e, Q')r\l + -^)-\ (1.7) 

and rewrite (1.6) as 

t" 4t ■ - (1 + + = 0. (1.8) 

Jo n + 1 (V n + 6 ^ 

where we have exchanged the order of the integrals in the second term of Eq.(1.6). 
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To determine the second unknown function /(^, Q^), we decompose it into two pieces 

f{^,Q') = f^'\^,Q') + f'\^,Q'), (1-9) 

where Q^) ~ and /(^(C, Q^) ~ 0{M^/Q^) is a small term. Since Q^) 

in the hmit — > oo, the function must satisfies large-Q^ behavior (1.5). We 

purposely choose 

Q') = (1 + + 26) ^i^. (1.10) 

From (1.8), (1.9), and (1. 10), we have 

4^ .e-i/<..K,<j=, . (1 ^ . 0. ,1.11) 

To determine small unknown function f^^\i, Q^), we put 

f'\^,Q')^^ [^^^^ dyrj(y,Q% (1.12) 
where r){y, Q"^) should be the order of 0{M'^ /Q^). Substituting (1.12) into (I-ll), one obtains 



r,«,<3^) = -3£tt,<3=)(l + ^)i^i^i^, (1.14) 



This equation can be satisfied for all n (n=2,4,...) only if the term in the bracket vanishes. 
Therefore 

which is indeed the order of 0{M^/Q'^). Substituting (1.9), (I.IO), (1.12), and (1.14) into 
(1.4), we finally obtain the modified W-W relation Eq.(4) with correction factors (5a), (5b), 
and (5c). 
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FIGURES 
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FIG. 1. The kinematical correction factors (see Eqs.(5a), (5b), and (5c)) plotted as functions 
of the Nachtmann variable ^ at Q^=3 (GeV/c)^. 
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FIG. 2. The correction ratios (see Eq.(6a,b)), plotted as functions ^, for Q^=3, 5, 10 and 
(GeV/c)2. 
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Solid line: gj^Cf), eq.(4) 
Dotted line: g2**(x), eq.(3) 
Dashed line: ga (f), eq.(6) 
Data at <Q^>=3.0 GeV*^ 
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FIG. 3. Modified #^(0, Eq.(4) and Eq.(6), and the original W-W relation, Eq.(3) plotted 
as functions ^, for Q^=3 (GeV/c)^. Data are taken from E143 and E155. 
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